New methods to determine the Hausdorff dimension of vortex loops in the 3DXY 

model 
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The geometric properties of critical fluctuations in the 3DXY model are analyzed. The 3DXY 
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I. INTRODUCTION 

It has been proposed that a second-order phase tran- 
sition is characterized by the breakdown of a generalized 
rigidity associated with a proliferation of defect struc- 
tures in the order parameter In the case of extreme 
type-II superconductors, it has been explicitly demon- 
strated that the system undergoes a continuous phase 
transition driven by a proliferation of closed loops of 
quantized vorticity (2. El- I n particular, it has been 
demonstrated that the connectivity of the vortex-tangle 
(the tangle of topological defects of the system) , changes 
dramatically at the critical temperature of the system 
In other words, it is possible to describe the su- 
perconducting phase transition in terms of the prolifer- 
ation of defect structures (viz. vortices and flux lines), 
which are singular phase fluctuations and determine the 
critical properties of the theory. Using this argument, it 
is possible to relate the critical properties of the phase 
transition to the geometric properties of the loops [5| at 
the critical point. 

The phenomenology of superconductivity usually 
starts with the Ginzburg Landau (GL) model. The 
Hamiltonian for this model is 



H{ W )=rnl\^ + ^-\^ + \D^ 



7^(Vx A) 2 .(l) 



Here ip = \tp\e l6 is a complex matter field, coupled 
to a massless gauge field A through the minimal cou- 



pling D u 



dn - iqA. 



is the mass parameter 



for the ip field, and is the self-coupling. This is a 
rich model which embodies many different aspects of 
super conductivity |y], in this paper we will focus on the 
limit k ^> 1 corresponding to extreme type-II supercon- 
ductors. In this case the gauge field fluctuations can be 



ignored, and we are left with a neutral theory. For 
this theory it can be shown that amplitude fluctuations 
in ip are innocuous 0, 3, and only the phase variables 
must be retained. When the resulting model is defined 
on a lattice, one obtains the 3DXY model. 

The 3DXY model is a model for phase fluctuations, and 
both the smooth spin wave fluctuations, and the singular 
transverse fluctuations are accounted for. The transverse 
phase fluctuations are defined by 



V x V0(r) = 27rn(r) 



(2) 
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where n(r) is the local vorticity. The vortices are the crit- 
ical fluctuations of the theory, which drive the superfluid 
density to zero in the neutral case, and the magnetic pen- 
etration length to infinity in the charged case. We have 
performed simulations directly on the phase degrees of 
freedom, and extracted the vortex content according to 
Eq. alternatively it is possible to integrate out the 
spin wave degrees of freedom and retain only the vortex 
degrees of freedom 7]. The definition Eq. © ensures 
that V x n = everywhere, hence the vorticity must be 
in the form of closed loops. 

Since the vortex loops are the critical fluctuations of 
the theory it is natural to formulate a theory expressed 
in terms of these degrees of freedom. In d — 3 it is pos- 
sible to start with the charged theory Eq. (|TJ in a fixed- 
amplitude approximation and deriveQ a field theory for 
the vortex loops. It turns out that this dual theory is a 
neutral \(j>\ 4 theory. Hence, the dual of a charged super- 
fluid is a neutral superfluid and vice versa, and further- 
more the vortex loops of a neutral (charged) superfluid 
are described by a field theory isomorphic to a charged 
(neutral) superfluid. In this paper we will use the conven- 
tion that ip represents the original superfluid, and that <p 
is the corresponding dual. Then (j> will be the field theory 
for the vortices of ip. 

The description in terms of loops is physically appeal- 
ing. Firstly, it highlights the physical meaning of the 
\(f)\ A term which, depending on the sign of 11$, represents 
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a steric repulsion (as in the case of type-II supercon- 
ductors) or a steric attraction for type-I superconductors 
with a first order transition^. In the case of a neutral 
superfluid the vortices will interact attractively through 
a long range interaction mediated by a dual gauge field, 
this will yield a vortex tangle more dense than a set of 
random loops 0- 

A geometric interpretation of the critical point in terms 
of proliferating geometric objects was given already in 
1967 by M. E. Fisher He considered droplets of one 
phase immersed in a background phase. Close to the 
critical point the distribution of droplet size was argued 
to behave as 



n(s) oc — e 
s T 



e oc IT - T n 



(3) 



n(s) is the mean number of droplets of "mass" s. The 
behavior of Eq. © is governed by two critical exponents 
a and r, where the former governs the vanishing tension 
when approaching the critical point and r is related to 
the entropy of a droplet. These two exponents can be 
related to the six ordinary critical exponents 



a = 2 



VG = 2 



T- 1 



d(T - 3) 

T - 1 



0G = S G = 1 + J 1 

a d{T — 2) 

3-r 



r - 1 
da 



7G = 



(4) 



The distribution n(s) Eq. J3J is also used to describe the 
cluster density close to the critical point in percolation, 
and the relations Eq. @ can be easily derived from that 
context as well 10] . The exponents in Eq. J2Jl have an 
index G to emphasize their geometric origin. These expo- 
nents will in general not agree with those of the underly- 
ing model. The Fortuin-Kasteleyn clusters ^lj of the Q- 
state Potts model is a special case where the exponents 
derived from r and a agree with those of the ordinary 
Potts model, i.e. (3q = (3, jq = 7, . . .. The vortices of the 
3DXY model, which we study in this paper, constitute 
another special case, here the exponents derived from r 
and a agree with the dual of the initial theory. Since the 
dual of a neutral superfluid is a charged superfluid, the 
study of the vortices in the 3DXY model can actually 
be used to glean knowledge of the critical properties of 
a charged superfluid. In particular, this can be used to 
relate the anomalous scaling dimension of the Ginzburg 
Landau theory, to the fractal dimension of the vortex 
loops in the 3DXY model. 

The paper is organized as follows. In section ^ we 
derive the relation between the fractal Hausdorff dimen- 
sion Dh of the loops and the anomalous dimension of 
the dual condensate, rj^. Sections IlIII and llVl are devoted 
to presenting our results from Monte Carlo simulations. 
As a benchmark, we determine the fractal dimension of 
percolation clusters in section ITTT1 and in section llVl wc 
determine the fractal dimension of the vortex loops of 
the 3DXY model. Finally, in Section [V] we discuss our 
results. 



II. RELATION BETWEEN THE HAUSDORFF 
AND THE ANOMALOUS DIMENSIONS 

The anomalous dimension r]^ of the <p held is defined as 
the critical exponent of the correlation function as follows 



G(x,y) = ^(x)0t( y )). 



(•5) 



This correlation function has the standard form at large 
distances at the critical point, namely 



G{x) ~ l/x d - 2+71 <» 



00, 



(6) 



where d the spatial dimension of the system. Particle- 
field duality dictates that this correlation function has a 
geometric interpretation, yielding the probability ampli- 
tude of finding any particle path connecting x and y. In 
the present work, the particle trajectories correspond to 
vortex loops which can be deduced from the phase distri- 
bution of the matter field. It is essential that any vortex 
path connecting two points x and y must be part of a 
closed vortex loop, since only closed vortex loops provide 
the vortex paths in this system. To highlight the impor- 
tant properties of the probability amplitude -P(x, y, N) 
of connecting two points x and y by a continuous vortex 
path, we start by focusing on the properties of the vortex 
loops. 

It is known that for random loops in a lattice, the 
number of steps N, defined as the number of occupied 

links, and the average distance Rp = y^A^ + + A§, 

where Aj are the coordinate variations of the trajectories, 
are related by 



R F ~ N 1/Dh 



N 



00 



(7) 



where Dh is the Hausdorff dimension . 
Since these vortex are the critical fluctuations of the the- 
ory their average size Rf is related to the correlation 
length £ of the field <fi so that we can set £ ~ Rf 



(8) 



With this we can now relate the correlation function 
G(x, y) with the Hausdorff dimension Dh by taking ad- 
vantage of the fact that in every second order transition 
the system is scale invariant at T = T c . We can, then, 
write the following scaling Ansatz for the probability am- 
plitude P(x, y, N) 



P(x,y,N) ex P(\r\,N) cx —^F 



H 



(9) 



where P(x,y,N) is the probability of finding a loop of 
length N connecting the points x and y. The probability 
of coming back to the starting point is generally given 
by P(x,x, N). For some models like self avoiding walks 
(SAW) and random walks in d > 3 this probability is 
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vanishing, and the scaling function F(z) has the limiting 
behavior 



lim ; 



F(z) 



i?> 0. 



(10) 



In the case of vortex loops, which are closed, we clearly 
must have i? = 0. This is discussed in more detail in 
Appendix ^ At the critical point the two-point corre- 
lation function G(r), Eq.lJ^jl, scales with an anomalous 
dimension rj^. At the same time we know that, by the 
definition of P(x,y,N) together with Eq.©: 



G(x, y)=J2 P & y> N ) = Yl N- d /° H F 



N 



N 



(.11) 



If we focus on the long loop regime (N 3> 1), we can 
replace the summation with an integral and obtain 



C 



\x-y\ 



l-D, 



-(12) 



Therefore, we can relate the anomalous dimension rj^ to 
the fractal dimension of the vortex loops, Dh, as follows 



V4> 



D 



H 



2. 



(13) 



We can also relate Dh and, indirectly 77^ , to another 
critical exponent, r, which is related to the density of 
loops in the system. From the definition of loop density 
we have 



D(JV)«l^P(x ia; ,iV) 



(14) 



here, P(x, x, N) is the probability of finding a closed loop 
of length TV and D(N) is the mean number of loops of 
length N. Near T c the loop distribution will then take 
the form 



D(N) ~ AT 



(15) 



If the system is homogeneous, all the contributions in 
Ea. (|14f> are equal. Thus, picking an arbitrary x we find 



N 

P{x,x,N) cx —D{N) oc N 1 ^. 



(16) 



Finally, combining Eq.@, Ea. (|16|) and Eq.{7|) we arrive 
at the following relations 



d 



1 



T-l 



(17) 
(18) 



It can be seen from Ea. (|13f, and Ea. l|18l) that the anoma- 
lous dimension of the dual condensate can be determined 
either through the determination of r or D^. the former 
being statistically easier to evaluate. However, 77^ is quite 
sensitive to r [drj^/dr — d/(r — l) 2 ~ 1.5]. 



Percolation Transition for 3D simple cubic SITE 
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Vol=40x40x40 
Vol=50x50x50 
Vol=60x60x60 



0.28 0.29 0.3 



conccntr LitLiin 



FIG. 1: (Color online) Percolation probability of sites as a 
function of concentration p for a simple 3D cubic system. The 
concentration at which percolation takes place is given by the 
common inflection point of the different curves for increasing 
lattice size. 



Percolation Transition for 3D simple cubic BOND 

Concentration threshold p =0.248627 
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FIG. 2: (Color online) Same as Fig0 but for a bond perco- 
lation 



The purpose of the present work is to carry out a direct 
evaluation of Dh which facilitates a check of the scaling 
law Eq. (|1 T|) . For this reason, the first part of this study 
has been devoted to test their validity and the validity of 
the algorithms used in the case of a percolation transi- 
tion, since it has been extensively studied in the literature 
and can be easily implemented. 

Before we start, let us return to the points made above 
concerning the HT-graph expansion of the partition func- 
tion and two-point correlation functions |13lll4| . Suppose 
one were to compare the Hausdorff dimension of physical 
vortex loops with the Hausdorff dimension of the sort 
of open ended, self-avoiding walks that one would con- 
sider in a HT-graph expansion of the two-point correla- 
tion function |13| ■ We may compute the loop-distribution 
function for vortex loops, D\ oop (N) ~ 7V~ Tlo °p, with 
the distribution of "closed" graphs in self-avoiding walks 
-Dsaw(^V) ~ A^~ rsAW . A self-avoiding walk of length N 
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is defined as "closed" if one after N steps ends up on 
a lattice point which is nearest neighbor to the starting 
point 0Q, see also Fig. 1 of [3. From Fig. 1 in Ref. 
[14j. it is immediately clear that the stop criterion for 
defining a "closed" self-avoiding walk, if applied to real 
physics vortex loops, would mean that one would discard 
a tail in the distribution function involving large vortex 
loops. This is so since one could connect the remaining 
last bridge by an arbitrarily long and complicated vortex 
path, not just over the shortest link. (In fact, these two 
point were connected by an arbitrarily long and compli- 
cated path, and there is no reason to close the path only 
by the shortest distance). This would lead to an overes- 
timate for Tjoop, since removing large loops means that 
D(N) drops too much as a function of N. An overesti- 
mate of r would lead to an underestimate of Dh , cf. Eq. 
p7[l. In the scaling relation rjj, + Dh — 2, it is crucial to 
compute Dh from the real physical vortex loops of the 
system. If one underestimates Dh by for instance con- 
sidering the fractal Hausdorff dimension of other objects 
than real physical vortex loops, this might lead one to 
erroneously conclude that one needs to introduce an ad- 
ditional exponent i? > and modify the scaling relation 
to rjrj, + Dh =2 — i?. While -d > is needed in the case 
of self-avoiding walks, this is not so for vortex loops due 
to the topological constraint V • n = 0. 



III. RESULTS FOR THE PERCOLATING 
SYSTEMS 

Percolation theory is used to describe a variety of natu- 
ral physical processes where disorder is an essential ingre- 
dient. Applications range from spontaneous magnetiza- 
tion of diluted ferromagnets, formation of polymer gels, 
electrical conductivity of amorphous semiconductors and 
many others [Tof. 

In link (site) percolation the links (sites) of a lattice are 
occupied independently with probability p, and then the 
resulting clusters of connected links (sites) are analyzed. 
For small p we will only have small clusters, and the 
there will be no path connecting the edges of the system. 
For large p there will be large clusters comparable to the 
entire system, and the system can sustain a current from 
edge to edge. The transition from an insulator to metal is 
a second order (geometric) phase transition at a critical 
value p c . 

Figs ^and show our results for the percolation prob- 
ability, i.e. the probability that there is at least one clus- 
ter spanning the whole system, for sites and bond per- 
colation respectively. We have found p c — 0.311619 and 
p c = 0.248627, respectively. These values agree reason- 
ably well with existing results Clusters have been 
identified using the Hoshen-Kopelman (HK) 0] labeling 
algorithm and up to 10 4 configurations were considered 
in the measurements. 

We next focus on the calculation of the fractal dimen- 
sion of percolation clusters at the critical point, and com- 



pare with results from the literature 17]. Formally the 
fractal dimension can be defined from the box counting 
technique |12|. The fractal object is covered with boxes 
of linear size and we count the number of boxes needed 
for complete coverage. The fractal dimension is inferred 
from the variation in total number of boxes N(l) with 
box size I as follows 

N(l) cx lim r° H . (19) 

Defining I on a lattice is difficult since the lattice spacing 
a is a fixed dimensionless constant. Instead, we can relate 
I to 1/Rp, which is equivalent to considering loops of 
different size Rf as the same loop, but at a different 
length scale. This relation is unique provided that the 
different loops can be embedded in a box A x A y A z of 
the same shape, i.e. 

^ = Cij V i + i i,j e [1,3]. (20) 

where Cy are constants enforcing a constant proportion 
constraint (CPC). This implies that loops have to be di- 
vided into groups, each one characterized by the actual 
values of the two independent ratios Cy, say C xy and 
C xz , and Ea. l|19|) applies separately to each group. In 
the current paper we have only considered cubic boxes 
with Aa; = Aj, = A z . 

From Eq.Q), it can be seen that there is a relation 
between the dimension of the cluster (i.e. the number 
of cells) and its radius (i.e. the dimension of the box 
containing it). Thus, we first relate the number N of 
occupied cells by the loops to their length Rp using an 
embedding box of fixed shape A^ = A y — A z . Although 
this approach is very close in spirit to the definition Eq. 
(|19fl the method fails miserably giving Dh = 1.329(8), 
which deviates strongly from the literature results Dh = 
2.5230[l7|. 

The problem, we believe, is that this method (later 
called Method!) is affected by two strong restrictions. 
Firstly, every box must scale in all the three directions in 
the same way (because of the CPC), i.e. all the boxes have 
to be cubes. Secondly, the extraction of the box can be 
mistaken if the cluster extends beyond the borders of the 
lattice, so that boxes touching the borders are excluded. 
For these two reasons the statistics of N(Rp) in Method! 
is extremely poor, and a correct value of Dh cannot be 
found. 

Another way to extract Dh is to consider the largest 
cluster for every configuration, and assume that the 
box containing this cluster consists of the entire lattice. 
When the linear size of the lattice is varied this gives the 
fractal dimension from 

N*(xL Dh , (21) 

where N* is the size of the largest cluster and L is the 
linear extent of the lattice. Using this approach we have 
found Dh = 2.52, which is consistent with the value 
t = 2.189. 
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FIG. 3: (Color online) Percolation probability Ol as function 
of temperature T for the 3DXY model at different lattice 
sizes. The critical temperature is given by the common in- 
flection point of the different curves for increasing lattice size. 
The value obtained is T c = 2.18727 

Hence, the study of the percolation transition has been 
useful to first test the simulation and the algorithms used 
for the extraction of the geometrical properties and, ul- 
timately, it has provided a test of Ea.l|17|l. We can 
now proceed to the more difficult task of determining 
the Hausdorff dimension for the ensemble of vortex loops 
generated by thermal transverse phase fluctuations of the 
3DXY model. 



IV. RESULTS FOR THE 3DXY MODEL 

The Hamiltonian for the case of a neutral condensate 
in the London approximation (i.e. neglecting amplitude 
fluctuations) is given by 

H(q = ) 6) = -J^2coa(e i -e j ). (22) 

Eas. l|17|l and (|13[l are valid only at the transition point 
where the system is scale invariant. Therefore the simu- 
lations must be restricted to the case T — T c . We have 
used the standard Metropolis algorithms for lattice sizes 
of L = 80, 100, 150, 180, all the runs had a hot start and 
5000 sweeps were discarded before measurements were 
made. Vortices are found identifying singularities in the 
phase configuration according to Eq. The temper- 

ature found (T = 2.18727, cf. FigOJl is close to, but 
different from the thermodynamical one, T c = 2.20184 
[l8|. Variations with the dimension of the system are 
still present, and it cannot be excluded that in the ther- 
modynamical limit the two values merge together. More 
extensive simulations along with a finite-size scaling anal- 
ysis will be necessary to provide this question with a more 
satisfactory answer. 



FIG. 4: Vortex-loop distribution D(N) ~ N~ T as a function 
of perimeter length N at the temperature T — 2.18727, as 
explained in the text. 



For the time being, since any T T c is better than 
T < T c to test Ea. l|17l) . all the following simulations will 
be done at the fixed temperature of T c = 2.20184. Now 
that the temperature has been set, r and Dh can be 
determined separately and their values compared to the 
ones available in the literature. 

All the simulations were performed on a system of fixed 
size L = 200 using up to 40 2Ghz Pentium4. The Monte 
Carlo runs had a cold start and 10 4 sweeps were discarded 
for thermalization before taking measurements over 4 - 10 5 
configurations. The loops in the system are identified 
with a modified version of the Hoshen-Kopclman algo- 
rithm which resolves vortex crossings randomly. The ex- 
ponent t for the loop distribution has been determined 
earlier as t — 2.310(3) 0,0. The results obtained for r 
are shown in Fig0J 

We now turn to the fractal dimension Dh- The first in- 
teresting fact is that the method used to extract the pair 
(length, box) used in the percolation transition does not 
work any more, giving unphysical values for Dh- On the 
other hand, the more correct method of extracting the 
exact box containing the loops yields a value in agree- 
ment with Ea. ljrTJl . 

There is also a topological reason for which we can 
expect to have different statistics for the two systems 
(and therefore we need to apply two different methods). 
When extracting the clusters in a percolating system, all 
connected sites (or bonds) are considered part of the same 
object, i.e. belonging to the same cluster. In general, this 
will lead to fewer but larger clusters in the system, this 
is particularly the case at the transition point, where the 
cluster tension vanishes. In the case of vortex-loops, on 
the other hand, two intersecting loops are connected with 
50% probability, which directly influences the value of r 
which is smaller in the case of percolation. 
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FIG. 5: Loop length N(Rf) as a function of the characteristic 
size Rf- The best fit has been performed only on the point 
between 4 < Rf < 47. 



FIG. 6: Variation of Dh with the lower cutoff. The satura- 
tion expected in the continuum limit is present but strongly 
effected by fluctuations. The finale value of Dh=2- 168(6) is 
obtained using a weighted fit over the last eight points 



Methodl 



Method2 



We have tested Methodl, and the extracted value of 
Dh = 2.168(6) (FigEJ is in good agreement with the 
expected one of Dh — 2.162(1), as calculated using 
Ea. (|17fl . and the value of r = 2.3876(8) previously ob- 
tained (Fig§J. 

Eq.JJJ) is assumed to hold only in the limit N — ► oo, 
so that small deviations from linearity are expected for 
small N where the discrete form of the lattice becomes 
important. Therefore, in order to reduce these lattice 
effects, we have varied the lower cutoff on Rp systemati- 
cally. Results are shown in Fig [13 which is characterized 
by a 'saturation' in the long-loop regime. To obtain a reli- 
able value of Dh this saturation region has to be reached, 
but, although this is the case, there are still considerable 
fluctuations. 

For this reason and owing to the results obtained for 
percolating systems, we have considered two new meth- 
ods in order to determine values of the Hausdorff dimen- 
sion Dh with better precision. The main limitation in 
Methodl is the limited statistics, in an attempt to im- 
prove on this situation we have devised two additional 
methods to calculate Dh- Both methods rely on 'relax- 
ation' of the CPC constraint. 

The strongest statistical limitation on Methodl is that 
the loops have to be divided into groups and that only 
comparison within each group are allowed. The conse- 
quence of this restriction can be directly seen in Fig|Sl 
where the longest loop considered has a linear dimension 
of ~ 47, to be compared with the linear dimension of the 
lattice of 200. Clearly then, the longest loops present in 
the lattice (the ones for which Eq. 0) holds) are excluded 
from the final statistics obtained. 



The first of these two new methods (Method2 below) 
determines a lower and an upper bound for Rp. To see 
how this can be done let us consider a box containing a 
loop with the condition = A y = A and A z = A + C 
where C is some constant, and A x , A„, A z the dimensions 
of the box containing the loop in the three directions. We 
immediately observe that, since C is constant, there is a 
violation of the CPC 



A 



A' 



A + C A' + C _ 

A, ~ A ^ A' ~ A 



(23) 



However, the error is of the order of C/A, so in the limit 
of large boxes (A 3> C) this approximation is exact. We 
can then define the following loop size Rp as follows 



R F = AV3 

R+ = (A + C)V3, 



(24) 

(25) 
(26) 



where R F and Rp are lower and upper bounds for Rf. 
By definition, we have Rp < Rp < Rp for every A, so 
that D H < D H < D H V A. Moreover, 



lim D~ = lim D+ = lim D+ = D H 

A— >OC A — >oo A— >oo 



(27) 



Then, D H and D H represent a lower and upper bond, 
respectively, for the Hausdorff dimension, thus allowing 
a direct estimate of the error in Dh- As can be seen from 
Fig0 the separation between the two values obtained, 
Dh- and Dh+, is still of the order of the fluctuations in 
Fig. hence there is no significant improvement in the 
Dh determination. 
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FIG. 7: Dependence of Dh- and Dh+ on the lower cutoff 
(Rf)- The deviation between the two values is still too large 
even in the long-loop regime. The finale value of Dh=2.14(7) 
is obtained averaging over the final values for D H and D H 



FIG. 9: Dependence of Dh on the lower cutoff (RF-min)- The 
saturation region is, as can be seen, much clearer compared 
with previous methods. The finale value of Dh ~ 2.1693(4) 
is obtained using a weighted fit over the last three points 
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R F 

FIG. 8: Loop length N(Rf) as a function of the character- 
istic size Rf- The best fit has been performed on the points 
between 4 < Rf < 150, to be compared with that of Methodl, 
in FigEl 



C. Method3 

We finally turn to the last method (Method3 below) 
which will prove to be the best one to extract a pre- 
cise value for Dh- In contrast to what is the case for 
Methodl and Method2, in Method3 there is no constraint 
on the proportions of the box containing the loops, since 
it is based on the approximation that, considering a large 
number of configuration, the different values for Dh will 
eventually compensate. As can be seen from FigEl this 
method gives much better statistics both in the small- 
and large-loop regimes. Moreover, the longer loops con- 
sidered in the statistic, having a linear dimension of 
~ 150 (cf. Fig|8}, although greater than the one con- 



sidered with Methodl and Method2, is still smaller then 
the linear dimension of the lattice. Hence, finite size ef- 
fects should not substantially influence the final value of 
Dh obtained. 



V. CONCLUSIONS 

The results obtained for the different methods are sum- 
marized in Tab[U 

All methods are in agreement with Ea. (|17fl and con- 
sistent with each other, although they have different ac- 
curacies. As predicted by ref.[19(, finite size effects are 
present both in the short and long- loops regime (Fig0J, 
but there is a clear monotonic behavior towards the con- 
tinuum limit (Figs EE and EJ. The best results are ob- 
tained from Method 3, which yields a final value for the 
anomalous dimension of the condensate?^ = —0.1693(4). 
This is consistent with both perturbative renormalization 
group calculations [2(3, El H2, and previous simula- 
tions HSU. 

In ordinary local quantum field theories without gauge 
fields, one can prove that r]^ must be greater or equal 
to zero. In contrast, in gauge theories the proof of a 
non-negative rj^ is not applicable due to the nongauge- 
invariant form of the correlation function [3 H3- El- 
Moreover, a negative value of r\§ implies a fractal dimen- 
sion larger than that of Brownian random walks, which 
means that the current lines are self-seeking. This point 
has been thoroughly discussed in Indeed a Dh > 2 
(corresponding to rj^ < 0) is necessary for the possibil- 
ity of the existence of a phase transition driven by a 
vortex-loop unbinding even in the presence of a finite 
external magnetic field 0, Q ■ This value is different 
from the one obtained for the 3DXY model, showing that 
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TABLE I: Values of r and Dh, it can be seen that Eg . 1171 
is satisfied for all methods, and that all the values of Dh ex- 
tracted are consistent with each other. The last two columns 
are a check of internal consistency. 





Dh 


T 


n+ = 2 - 


r/4, = 2 — Dh 


Methodl 


2.168(6) 


2.3876(8) 


-0.162(1) 


-0.168(6) 


Method2 


2.14(7) 


2.3876(8) 


-0.162(1) 


-0.14(7) 


Method3 


2.1693(4) 


2.3876(8) 


-0.162(1) 


-0.1693(4) 



the gauge fluctuations of the Abelian Ginzburg-Landau 
model (dual 3DXY) modify the critical behavior, so that 
the Abelian gauge model and the </> 4 model belong to two 
different universality classes. 
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APPENDIX A: VORTEX LOOPS AND EQ. ([Toll 

The quantity 

lim F ( O ) = Urn F(z) = z* > (Al) 

is a scaling function for the probability of traversing the 
distance between two points along a continuous vortex 
path, and holds for an ensemble of both closed and open 
self-avoiding walks |27j . Eq. HAlj) is often used in polymer 
physics to describe the statistical properties of polymer 
tangles. The statement that d > is certainly not 
true for the statistical properties of vortex loops, which 
are, by construction, closed and not self avoiding, so we 
have P(x,x,N) ^ and ■& = [28|. The main physical 
difference between vortex loops in a superconductor and 
a tangle of polymers, is that while polymers can start 
and end inside a system, this is not so for vortex loops. 
They must either close on themselves or alternatively 
thread the entire system. If we denote by n the local 
vorticity arising from transverse phase fluctuations, 
defined by Eq. J2J , then we have the constraint V • n = 



everywhere inside the superfluid (and superconductor). 
We reemphasize that this in particular means that 
each and every vortex path in the system must be part 
of a closed vortex loop. No paths of vortices can be 
open-ended. Hence, for the former we can have d > 
while for the latter i? — 0. To emphasize this further, 
let us consider the physical meaning of a positive •& 
[l4j . In the problem of self-avoiding walks, ■& > is 
an exponent that governs the asymptotic number of 
open-ended paths of length N at the critical point, via 
the relation zn ~ N^/ Dh . We also emphasize that 
in this paper and in previous works, we have 0, 0, Q 
exclusively been dealing with the geometrical properties 
of vortex-loop paths at the critical point in supcrfluids 
and superconductors. 

Let us also make another remark which is relevant in 
this context. It is known p9j. that the partition func- 
tion of the 3DXY model may be expanded in so-called 
high-temperature (HT) graphs, which are isomorphic to 
the vortex-loop gas of a superconductor. In the case 
of charged superconductor the fluctuating gauge field 
screens the Coulomb interactions, and the resulting loop 
interactions are short range Q. Hence, for the purposes 
of studying the critical thermodynamics of a supercon- 
ductor, one might as well use such a HT expansion of 
the partition function, instead of summing over all con- 
figurations of closed vortex loops of the system. In other 
words, at the level of the partition function, one can ele- 
vate the HT graphs to objects that are equivalent to real 
physical topological defects in the superconducting order 
parameter, i.e. closed vortex loops. A similar property 
is known for the two-dimensional Ising model, where one 
at the level of the partition function can proceed with a 
HT expansion and identify a formal equivalence between 
the closed HT graphs and the topological defects of the 
theory, namely closed lines in a 2D plane connecting do- 
mains of oppositely directed Ising spins. 

However, when we compute correlation functions, as 
we need to do for computing the probabilities we have 
discussed above, the connection between physical closed 
vortex loop paths and the graphs of the HT expansion 
for the correlation functions, is far less obvious. When 
one utilizes the HT expansion to compute the correla- 
tion function ((/>(x)(^(y)), one has to consider open-ended 
graphs starting at x and ending at y |14j. One may, if 
one wishes, define a closed such path as the path at which 
the endpoint of the walk for the first time reaches a lat- 
tice point which is the nearest neighbor of the starting 
point, see for instance Fig. 1 of |14| . One may further 
go on to discuss the geometrical properties of such open- 
ended (and "closed") HT graphs involved in computing 
the two-point correlation function in a HT expansion. 
However, the graphs that one then ends up with study- 
ing, have nothing to with the closed vortex paths that 
are the topological defects of the superconductor, and 
which were the objects under consideration in our pre- 
vious work [5j. (The fractal structure of HT graphs in 
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O(N) models in two spatial dimensions, as well as the 
fractal structure of spin clusters and domain walls in the 
two-dimensional Ising model, has recently been investi- 
gated in detail [H [Si.) The Hausdorff dimension of 
suchcorrelation function HT graphs will not be the same 



as the Hausdorff dimension of the closed vortex loops of a 
superconductor. Such open-ended graphs violate the con- 
straint that the vortex-loop paths must respect, namely 
that they must form continuous closed paths and cannot 
start or end inside the system. 
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